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We investigate in detail the density dependence of the symmetry 
energy in a relativistic description by decomposing the iso-vector 
mean field into contributions with different Lorentz covariant prop- 
erties. We find important effects of the iso-vector, scalar channel 
(i.e. (5-meson like) on the high density behavior of the symmetry 
energy. Applications to static properties of finite nuclei and to dy- 
namic situations of heavy ion collisions are explored and related 
to each other. The nuclear structure studies show only moder- 
ate effects originating from the virtual 5 meson. At variance, in 
heavy ion collisions one finds important contributions on the reac- 
tion dynamics arising from the different Lorentz structure of the 
high density symmetry energy when a scalar iso-vector 5 field is 
introduced. Particularly interesting is the related neutron/proton 
effective mass splitting for nucleon transport effects and for reso- 
nance and particle production around the threshold. We show that 
the (5-like channel turns out to be essential for the production of 
pions, when comparing with experimental data, in particular for 
high momentum selections. 
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1 Introduction 



Heavy ion collisions at relativistic energies from 0.1 up to 1 — 2 AGeV offer 
the possibility to access the equation of state (EOS) of nuclear matter under 
extreme conditions of density, isospin and temperature [1]. This knowledge is 
important in understanding many interesting astrophysical phenomena such 
as the physical mechanism of supernovae explosions and the neutron star 
structure and cooling. 

During the last three decades one has attempted in many ways to investigate 
the properties of highly excited hadronic matter. In phenomenological models 
of a relativistic quantum field theory, like the Relativistic Mean Field (RMF) 
model of Quantum-Hadro-Dynamics (QHD)[2], a significant decrease of the 
effective nucleon mass with density was found to be essential for a reliable de- 
scription of infinite nuclear matter and finite nuclei. More microscopic studies, 
the so-called Dirac-Brueckner-Hartree-Fock (DBHF) approach, based on real- 
istic nucleon-nucleon (NN) interactions in the One-Boson-Exchange (OBE) 
model, see refs.[3-6], also predict a strong decrease of the effective nucleon 
mass. Similar effects have been found for particles with strangeness such as 
kaons [7-9] . The results can be represented by saying that the effective nucleon- 
meson vertices or coupling functions exhibit a density dependence and, for en- 
ergies higher than the Fermi energy, an additional momentum dependence [6] . 
As a consequence the EOS predictions of different models show quite variable 
density dependences for densities higher than saturation, characterized by re- 
ferring to a hard or soft EOS, which have been tested in heavy ion collisions 
[10]. 

So far asymmetric nuclear matter has been only poorly investigated, but cer- 
tainly the same considerations appear for the density dependence of the sym- 
metry term which are directly related to the Lorentz structure of the effective 
fields in the isovector channel and to their weights (the coupling functions) in 
different baryon density regions. Around the saturation point the 0,4 parameter 
of the Weizsaecker mass formula contains the leading information (in a sense 
equivalent to the a± parameter for the isoscalar part). Actually, there is quite a 
wide range of extracted values, roughly from 28 to 36MeV, depending on the 
way the surface effects are accounted for. This clearly shows that even to re- 
produce the equilibrium EOS point we need a better knowledge of the density 
dependence of the symmetry term since the slope (the symmetry pressure) 
and the curvature (the symmetry compressibility) around p are also playing 
an important role in reproducing nuclear masses. A related sensitive observ- 
able appears to be the difference between neutron and proton r.m.s. radii in 
charge asymmetric nuclei. The wide range of predictions made by different 
relativistic (and non-relativistic Skyrme-if-F) effective interactions has been 
nicely presented recently in the refs. [11,12]. Other independent information on 
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the symmetry energy in dilute matter, below saturation, have been suggested 
from heavy ion collisions in the Fermi energy domain, in particular from the 
isospin distillation mechanism in fragment production, [13] and refs. therein, 
and from the isospin content of fast particle emissions, [14]. 

In standard RMF models the iso- vector part of the mean field is described in 
terms of a constant vector, isovector, p-coupling giving an almost linear den- 
sity dependence of the symmetry energy. Recently the importance of a con- 
tribution of the scalar isovector channel, the virtual 5[a (980)] meson, within 
the same RMF framework, has been stressed [15,16]. This mostly affects the 
high density behavior, as we will show in the following. We note that the 
5— channel usually has not been considered before just on the basis of the 
weak contribution to the free Nucleon-Nucleon interaction. But in the spirit 
of the Effective Field Theory as a relativistic Density Functional Theory, 
(the EFT/DFT framework, see [17,18]), the relevance of this channel could be 
completely different in nuclear matter, due to medium and many-body effects. 
We like to note that very recently , see the conclusions of ref.[12], the ^-channel 
has been reconsidered as an interesting improvement of covariant approaches, 
in the framework of the EFT/DFT philosophy. One of the main tasks of our 
work is just to try to select the dynamical observables more sensitive to it. 

In the Density Dependent Hadronic (DDH) field theory [19], the detailed den- 
sity dependence of the nucleon-meson couplings has been taken either directly 
from DBHF calculations [20-22] or has been fitted to data of nuclear matter 
and finite nuclei [23]. Certainly the correct microscopic approach should be 
to derive the coupling functions, in a QHD mean field scheme, from Dirac- 
Brueckner-Hartree-Fock calculations. Several attempts have been performed 
so far, see refs. [20-22], but the results do not agree and appear to be model 
dependent, in particular for the isovector channel. We will discuss such point 
in more detail below. 

Aim of the present work is to try to shed some light on this intricate puzzle 
of the baryon density dependence of the symmetry term, in particular in the 
high density region. For densities beyond saturation one needs extrapolations 
which can only be investigated in dynamical situations of heavy ion collisions 
at intermediate energies, where hot and dense asymmetric nuclear matter is 
formed. Indeed, the study of the high density symmetry energy is an object of 
recently started investigations [24-27]. Their motivation has been the planning 
of new experimental heavy ion facilities with neutron-rich radioactive beams 
with data available in the near future. 

Here we will actually explore the properties of asymmetric nuclear matter 
in different nuclear systems, i.e. finite nuclei and heavy ion collisions. For 
our studies we apply various versions of QHD models, in the Relativistic 
Mean Field (RMF) approach, which differ in the treatment of the effec- 
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tive meson-nucleon vertices. In the standard RMF models the nucleon-meson 
couplings are given by constants fitted to nuclear matter saturation proper- 
ties. Some non-linear contributions for the meson fields are considered, NL — 
RMF approaches, which imply definite density dependences in the model. The 
RMF — DDH theory uses the so-called Vector Density Dependent (VDD) 
parametrizations [19,23] in terms of directly density dependent effective cou- 
pling functions. In this case there is no need of non-linear terms. 

In this context relativity is important because it naturally explains essential 
nuclear matter properties such as the cancellation of large Lorentz vector and 
scalar fields, in the isoscalar channel, leading to a small value of 16 MeV 
for the binding energy per nucleon and, simultaneously, to a strong spin-orbit 
potential. Within a covariant formulation one is also able to better understand 
the origin of the repulsive and attractive contributions to the effective NN 
interaction. The separation of the total nuclear mean field potential in terms of 
its Lorentz components is of particular importance when studying dynamical 
situations of heavy ion collisions at relativistic energies [25,28,29]. 

In the RMF picture the nuclear mean field is given by the exchange of mesons 
with different Lorentz properties. The iso-scalar part which describes the EOS 
of symmetric nuclear matter is mainly characterized by an attractive scalar 
field (a-meson) and repulsive vector field (cu^-meson) both with contributions 
of the same order 400 — 500 MeV to the total energy per nucleon. The corre- 
sponding masses, fixed as m a ~ 550 MeV and m w = 783 MeV, are important 
in finite nuclei calculations. One should note the importance of both these 
mesons in describing the saturation properties of symmetric nuclear matter. 

A covariant treatment of the problem is also helpful for the understanding of 
the properties of asymmetric nuclear and neutron matter. As in the iso-scalar 
case, the iso-vector part of the nuclear mean field, which is responsible for the 
density dependence of the symmetry energy, can again be characterized by 
iso-vector mesons with different Lorentz properties. These are the iso-vector, 
vector field (p^-meson) and the iso-vector, scalar field (5-meson), with masses 
m p = 769 MeV (763 MeV in VDD models) and m s = 980 MeV. Due to less 
available experimental information several possibilities have been proposed to 
describe the density dependence of the symmetry energy. 

In contrast to the iso-scalar case, since here we do not have the stringent 
saturation condition of balancing attractive (scalar) and repulsive (vector) 
contributions, one can apply either only the p-meson, or both, the p- and (Ja- 
mesons. Because of the different Lorentz properties of the mesons one expects 
different density behavior of the symmetry energy especially at high densities. 
The argument is very simple and clearly discussed in the refs. [15,16]. Vector 
and scalar effective fields will contribute to the symmetry term with quantities 
proportional to baryon and scalar densities, respectively. As a consequence the 
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main differences should appear at baryon densities above saturation. Of course 
such direct interplay between p and 5 contributions to the symmetry energy 
can be reproduced in the symmetry term just by introducing density depen- 
dent iso vector /vector couplings without the need of a new virtual (scalar) 
meson in the isospin channel. The scalar nature of the 5-meson, however, also 
leads to a new effect, the splitting of the neutron and proton effective masses, 
with interesting consequences on nucleon transport properties and on particle 
and resonance production in the reaction dynamics at high baryon and isospin 
densities. In this work we will focus our attention on these new features of a 
isovector /scalar coupling in realistic collision simulations and on the related 
most sensitive experimental observables. Charged pion yields appear to be a 
suitable quantity to measure in order to see such effects. 

As we shall show, finite nuclei studies cannot separate easily between the 
effects originating from the different iso-vector mesons for densities mainly 
below saturation, see also [23]. Therefore, relativistic heavy ion collisions rep- 
resent a unique tool to study the properties of highly compressed asymmetric 
nuclear matter using radioactive beams. 

The paper is organized as follows. In Section 2 we discuss asymmetric nuclear 
matter within the RMF theory. The models are then applied to static cases 
of finite nuclei (Sec. 3) and to dynamical ones of heavy ion collisions (Sees. 4, 
5). The application to finite nuclei is important to demonstrate the reliability 
of the models in this sector, and the role of the 5-meson will be particularly 
emphazised. Then we describe briefly the theoretical treatment and numer- 
ical procedure of heavy ion collisions. The transport theoretical analysis is 
mainly applied here to particle production. In particular we focus the atten- 
tion on charged pion yields, with rapidity and transverse momentum selections 
(Sec. 5). A comparison with experiments is performed, when data are available, 
also in a preliminary status. Final remarks and an outlook follow in the last 
Section. 



2 Equation of state of asymmetric nuclear matter 

The starting point of a covariant quantum field theory is the Lagrangian of an 
interacting many body system of baryons and mesons. The latter determines 
the interaction between baryons in terms of nucleon-meson coupling vertices. 
Within a Non-Linear Relativistic Mean Field (NL — RMF) approximation the 
baryons are given by Dirac spinors and the mesons are described by classical 
field equations 

[-y^d" - &m,7° " 9 P l°T3Po - (M - g a a - g S T 3 6 3 )]V = (1) 
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15.60 


10.50 


1.10 


0.0 


0.015 


-0.004 




NL P 5 


15.60 


10.50 


3.15 


2.4 


0.015 


-0.004 




NL3 


15.73 


10.53 


1.34 


0.0 


0.01 


-0.003 





Table 1 



Coupling parameters in terms of = (£]-) for i = a, u, p, S, A = and B = 
for the non-linear RMF models using the p and both, the p and 5 mesons for the 
characterization of the isovector mean field in comparison with the iVL3 model. See 
text. 

m 2 a a + aa 2 + ba 3 = g a <tv >= g a p s (2) 
ml^ = g w < f >= g u f (3) 
™ 2 P P = 9 P < t 7 r 3 ^ >= 9 P P3 (4) 
m 2 s S 3 = 9s < *r 3 * >= g s p s3 ■ (5) 

The nuclear system is characterized by a baryonic quantum field and a sur- 
rounding background of classical mesonic fields, the iso-scalar, scalar a-, the 
iso-scalar, vector u>-, the iso-vector, scalar 5- and the iso-vector, vector p- 
mesons. The isospin vector and scalar densities are given by p 3 = p p — p n , p s3 = 
Psp — Psn, respectively, p Pin being the proton and neutron densities. The iso- 
scalar a field contains non-linear contributions. In DDH — RMF approaches 
the above equations will be modified for the presence of rearrangement terms 
in the nucleon self-energies and for the lack on non-linear terms, see refs. 
[19,23]. 

The energy-momentum tensor within the same non-linear RMF approach is 
given by (f { = (gi/rrii) 2 ) : 



+ i<T [m^ 2 + 2U(a) + f s p 2 s3 - f^j a j a - f pJ %j 3o 



(6) 



with the kinetic part, t^ v (x) = i < ^^d v i> >, and U(a) the non-linear 
contributions of the iso-scalar cr-meson U = |a<r 3 + \ba A . The effective mass is 
now different for protons and neutrons due to the appearance of the iso-vector, 
scalar 5-meson [30]: 



m* = M — g a a ± g$5 3 (- proton, + neutron) 



(7) 



For the investigation of asymmetric nuclear matter and neutron matter we 
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introduce the asymmetry parameter a = p "+ p p p which describes the relative 
ratio of the neutron to proton fraction of the nuclear matter. The asymmet- 
ric nuclear matter calculations can be done by solving selfconsistently the 
Eqs. (1-5,7). The parameters of the model, fixed to nuclear matter saturation 
properties, are given in the Table 1 and compared to the widely used NL3 
parametrization [31,32]. In order to isolate the effects in the isovector channel 
we have chosen on purpose very similar parameters in the isoscalar part. At 
saturation in symmetric matter we have a density p sa t = 0.148/m -3 , a nucleon 
effective mass m* N = 0.6M and a compressibility K NM = 211.1MeV . In NL3 
the symmetry parameter is a 4 = 37 AMeV while in our NLp, NLpS choices 
we fit a value a 4 = 30.5MeV. 

In all the above NL — RMF schemes the nucleon-meson couplings are given by 
constant values, see Table 1. We also want to compare to RMF models based 
on the Density Dependent Hadronic (DDH) field approach [19], since they 
have been shown to reproduce well nuclear matter and finite nuclei [20,23,33]. 
We will follow here in particular the parametrization of ref . [23] , where the den- 
sity dependence has been fitted to finite nuclei, and which makes a prediction 
for the isovector channel at high density. The symmetry energy at saturation 
is chosen as a 4 ~ 33 AMeV , see also ref. [33]. The corresponding properties 
for symmetric matter are: density p sat = 0.153/m~ 3 , nucleon effective mass 
m* N = 0.55M and compressibility Knm = 2A0.0MeV . So the interactions with 
density dependent couplings used in this work are softer with respect to the 
NL — RMF models. We will see that this will affect the particle production 
yields. 

However, the work of ref. [23] made no predictions for the 5-like field, since 
there was not much sensitivity to it in finite nuclei, as expected (see Sect.l) 
and as we will also see here later in Sect. 3. To test this part in heavy ion 
collisions we use three variants of these DDH models. DDHp is just the 
VDD parametrization of ref. [23], with one a p-like isovector field with an 
exponential density dependence of the coupling. In the models DDH3 we use 
the functional form of the vertex density dependences for the p- and 5-fields 
as extracted from the DBHF calculations of ref. [34]. In DDH3p we use only 
the p-field and normalize to the a 4 ~ 33AMeV value at saturation, as in the 
DDHp (VDD) case. In the DDHp5 we use p- and 5-fields and normalize, see 
following, as to obtain the same symmetry energy at saturation. 

The density dependence of the couplings, for all channels, for all the tested 
models is shown in Fig.l. In the NL — RMF models the couplings are con- 
stant, for the DDH — RMF model they are chosen as discussed above. Note, 
in particular, the unusual density dependence of the 5-like vertex and the at- 
tempt to reproduce some average constant value in the corresponding NLpS 
parametrization. Finally we like to remind that, due to the non-linear self in- 
teraction of the cr-meson, in the NL — RMF models the effective g a will also 
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show a density dependence. 



The symmetry energy E sym is defined from the expansion of the energy per 
nucleon E(p B ,a) in terms of the asymmetry parameter (we use p B for the 
baryon density in order to avoid confusion with the p meson): 

E(p B , a) = E(p B ) + E sym {p B )a 2 + 0(a 4 ) + ■■■ (8) 



with the abbreviation 



l d 2 E(p B ,a) I j?e 

*svm- 2 da 2 l«=o- 2 pB dpl lp3=0 ' [ ) 



From the definition of the energy momentum tensor (6) one obtains for the 
symmetry energy [15,16]: 



E = - — 



J sym 



6£ F 2 



fp-f'is;) 



PB 



(10) 



where -Ep = \jk F + m* 2 . 

It is seen that including the scalar iso-vector 5-meson the bulk asymmetry pa- 
rameter is given by the combination [f p — f 5 (^) 2 ] of the repulsive vector (p) 
and attractive (5) iso-vector couplings. We thus have to increase the p-meson 
coupling when including the 5-meson in the iso-vector part of the equation 
of state in order to reproduce the same bulk asymmetry parameter a^. This 
is a very general feature, not depending on the model, that we will see to be 
important in the dynamical simulations. Another 5-effect, of interest for trans- 
port properties, will be the splitting of the neutron/proton effective masses 
[15,16], see Eq.(7). 

From Eq.(10) we note the moderate role of the scalar, iso-vector 5-meson 
at low densities, and therefore on properties of finite nuclei (see below), as 
already pointed out in [23]. It is important to realize the different Lorentz 
structure of the iso-vector mesons. The vector p-meson contribution depends 
linearly on the baryon density whereas the scalar 5-meson part is proportional 
to the scalar density or correspondly to the factor (m* /Ep), decreasing at 
high baryon density. 

In Fig. 2 the density dependence of the symmetry energy, Eq. (10), for the 
different models used here is shown. The kinetic and the potential contribu- 
tions to the total symmetry energy can be also seen separately on the bottom. 
Within the Walecka models (NLp and NLp8)) one observes a similar behavior 
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Fig. 1. Density dependence of the different coupling functions as indicated for the 
relativistic mean field models used here (see text for details). 

of E sym at densities below saturation. With increasing baryon density ps, how- 
ever, the differences arising from the presence of the 5-meson in the iso-vector 
channel become more pronounced, as expected. This is due to the saturation 
of the scalar density (the factor (^) 2 of Eq.(10)) for high densities and, on the 
other hand, to the enhanced value of the vector p-meson contribution (needed 
to reproduce the same 04, see Table 1) which then grows linearly with baryon 
density. 

In the DDH models the high density influence of the 5-meson on the iso- 
vector part of the equation of state is more involved. We have now a rapid 
decrease with density of the p— coupling function while the increase of the 
5— coupling above saturation almost counteracts the relativistic scalar density 
effect, see Fig.l. As we can see from the two DDH3 columns of Fig. 2, with the 
inclusion of the 5-meson the symmetry energy is not much changed for baryon 
densities higher than saturation and actually can become even softer. The 
Lorentz structure of the isovector contributions is now different, and moreover 
we have an important splitting of the n/p effective masses. We will show that 
from particle production it may be possible to trace back these fine relativistic 
effects, although of course difficult. The density behavior of symmetry energy 
for the the DDHp model, (last column of Fig. 2), can be regarded as very soft 
originating from a relatively strong density decrease of the p-meson coupling, 
see Fig.l. This is particularly evident in the potential contribution (dashed 
line on the bottom-right). From Fig. 2 the importance of relativistic features 
in the understanding of the high density behavior of the symmetry energy is 
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Fig. 2. Density dependence of the symmetry energy, Eq. (10), for the same models as 
in Fig. 1. The lower part shows separately the kinetic (solid) and potential (dashed) 
contributions to the total symmetry energy. 

evident. 



In contrast to the iso-scalar part of the mean field, which has been extensively 
investigated in the last three decades theoretically and experiment aly [1,35], 
there is not enough empirical information about the high density behavior of 
the symmetry energy. Heavy ion collision experiments with radioactive beam 
facilities have recently been started or are being planned. In the following 
section we investigate first in finite nuclei the Lorentz structure of the iso- 
vector channel in the framework of the models discussed here. In the next 
section we will then pass to heavy ion collisions at intermediate energies, where 
higher densities can be reached and relativistic effects can be enhanced. We 
will discuss constant coupling Walecka-type models (NLp and NLpS)) as a 
kind of reference to clearly see the physics of the relativistic effects of different 
isovector contributions. We will compare then with the more sophisticated 
DDHp and DDH3p, DDH3p5 models in order make predictions for the main 
observables in experiments. 



3 Finite nuclei (static case) 



In this Section we focus on the role of the iso-vector, scalar 5-meson in the 
properties of finite nuclei and, particularly, we compare with other models 
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Fig. 3. Mass dependence of the binding energy per nucleon for Ca, Ni, Sn and Pb 
isotopes for the different models. 

widely used in nuclear structure studies. For the description of finite nuclei we 
solve the equations of motion for nucleons and mesons, Eqs. (1-5), including 
the gradient terms in the classical meson field equations in a selfconsistent 
scheme. In case of density dependent coupling functions also additional rear- 
rangement contributions of the iso-scalar and iso-vector classical fields (not 
shown in Eqs. (1-6)) are explicitely taken into account [23,19]. 

Fig. 3 displays the mass dependence of the binding energy per nucleon for 
different isotopes, starting from the most symmetric case (Ca) to the most 
asymmetric one (Pb). The results for all models are close together but the 
differences originating from the 5-meson become more pronounced with in- 
creasing neutron to proton asymmetry. This is particularly evident for the 
constant coupling RMP models, when comparing NLp and NLpS results. 

We recall from Fig. 2 that in the region around and below saturation, tested 
in finite nuclei, no significant differences between the models are seen in the 
symmetry energy. Only between the two extreme cases of NLpS and DDHp 
there is a visible difference in the low density regime which also leads to 
effects in the binding energy per nucleon for highly asymmetric finite systems. 
Generally, the binding energy becomes larger with the 5-meson included in 
the iso-vector mean field because here the attractive mean field is stronger 
(smaller) for neutrons (protons), with an observed net effect for neutron-rich 
nuclei. The comparison with the widely used models DDHp [23] and NL3 
[31,32] (here essentially given by the NLp points), which nicely describe finite 
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Fig. 4. Mass dependence of charge proton radii for Ca— , Ni—, Sn— and 
Pb— isotopes: (grey squares) exp. data, (full triangles) DDHp, (open triangles) NLS, 
(full circles) NLp) and (full squares) NLp5). 



nuclei properties, indicates the consistency of our parametrizations, that we 
will then extend to the description of heavy ion collisions. 

Fig. 4 shows the charge proton radii for different isotopes, as indicated. First, 
one observes moderate differences between the different models, since the pro- 
ton density profiles are very similar, as expected. All the models can reproduce 
the experimental values, as an important issue. 

One should expect larger effects of the internal structure of the symmetry 
energy by looking at the difference between neutron and proton radii, which 
is more sensitive to slope and curvature of the symmetry energy around sat- 
uration [11,12]. This is shown in Fig. 5 for the same models and systematics. 
As in the previous figures, the influence of the 5-meson is again not much 
pronounced in the NL— and DDH— models, because of the similar density 
behavior of the symmetry energy at low densities. We can see a systematic 
smaller r n — r p for neutron-rich nuclei, i.e. a reduced neutron diffusivity, in 
the case of DDH parametrizations, due to the larger p- (symmetry energy) 
repulsion seen by neutrons on the surface (at lower densities), see Fig.l. This 
seems to be in better agreement with data, see [12]. We note that this result is 
present also when the 5 channel is added, see the open squares of the DDH3p5 
interaction, in particular for the largely discussed 208 Pb case. 

At variance the slightly smaller symmetry repulsion for neutrons at densities 



12 





NLp 
m- -■ NLp5 

A- ■▲ DDHp 
O-O DDH3p 
H3p5 




170 180 190 200 210 220 230 
A 



Fig. 5. Mass dependence of the difference of the neutron to proton radii for the same 
isotopes and models as in Fig. 3. 



below p in the NLpS parametrization can consistently account for the larger 
neutron diffusion and the overall larger binding energies for very neutron-rich 
nuclei, Fig. 3. 



This discussion confirms the relatively small effect of a more extended rela- 
tivistic structure of the symmetry term in the low density range, tested in 
finite nuclei structure calculations, see also [23]. Similar results have been ob- 
served in the spin-orbit splitting of single particles levels, which are not shown 
here. 



Generally the properties of finite nuclei can be well reproduced in the NL 
models and are close to the results of more realistic DDH approaches [23]. 
This is an important test of the RMF parametrizations before applying them 
to the dynamics of heavy ion collisions. However, it seems difficult to determine 
the relativistic structure of the iso- vector part of the nuclear mean field in the 
static case, because the influence of the 5-meson is only moderate, except if 
one goes to extremely asymmetric neutron-rich nuclei. However, this is not the 
aim of the paper, but the behavior of the symmetry energy at supra-normal 
densities. Since in heavy ion collisions highly compressed matter is formed, at 
least for short time scales, we study now the dynamical case in more detail. 
The high velocity fields will also enhance the different Lorentz properties of 
the various contributions, [25]. 
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4 Heavy ion collisions (dynamical case) 



The theoretical treatment of heavy ion collisions at relativistic energies be- 
tween 0.1 and 1 — 2 AGeV is based here on a covariant transport theory. As 
detailed described in Ref . [36] , starting from a covariant quantum description 
of the nuclear many body problem one arrives at a transport equation which 
is called in the literature Relativistic-Boltzmann-Uhlenbeck-Uehling (RBUU) 
equation 



with effective masses to*, effective momenta k*^ and the field tensor defined 
by = - d^YT. are the vector self-energies and W(kk 2 \k 3 k 4 ) the 
transition rate for the process k + k 2 — > k 3 + k±. 

The main approximations involved in Eq. (11) are a semi-classical gradient ex- 
pansion and a quasi-particle approximation which sets the effective momenta 
on shell, i.e. A;* M A;* = m* 2 . The numerical procedure is performed within the 
relativistic Landau Vlasov method [37,38], where the phase space distribu- 
tion function f(x, k*) is characterized by test particles of a covariant Gaus- 
sian shape in coordinate and momentum space. Sufficient stability is achieved 
when using ~ 40 — 50 test particles per nucleon for the simulations. 

Protons and neutrons are propagating separately according to their hadronic 
fields and the Coulomb interaction. The collision integral is treated by Monte- 
Carlo methods [39] including all important inelastic channels up to 1—2 AGeV 
[40]. Here we explicitely propagate A 0,± ' ++ — and iV*-resonances including 
their decay to one- and two-pion (n ^) channels. At these energies we consider 
the lower mass resonances only, with free mass assignements Ma = 1232 MeV 
and Mat. = 1440 MeV. For elastic scattering we apply the empirical evalu- 
ations of isospin dependent cross sections [41]. For the inelastic channels we 
use the free cross section parametrizations of ref. [42]. The produced pions are 
also propagated under the influence of the Coulomb field. We should mention 
that the parameters in the collision term are commonly used by other groups. 
Main properties as rapidity distributions and transverse momentum spectra 
of protons as well as pion multiplicities are in agreement with results from 
other groups when using similar parameters for the nuclear equation of state 



k*»d x + {k* v F^ + m*d%m*) d k * f{x, k*) = 




-W(kk 2 \k 3 h)(2ir) 4 5 A (k + k 2 -k 3 - h) 



x " f{x, k 3 )fix, k 4 ) (1 - fix, k)) (1 - fix, k 2 )) 
fix, k)fix, k 2 ) (1 - fix, k 3 )) (1 - f(x, k 4 )) 



(11) 
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Fig. 6. Time dependence of the (-)-ratio (left), the multiplicities of A ,_ , A +,++ 
resonances (middle) and the (7r _ /7r + )-ratio (right) for central (b < 2 fm) Au + Au 
collisions at 0.6 (top) and 1 AGeV (bottom). Results of RBUU calculations with 
the p-meson (NLp, solid) only and both, the p and 5 mesons (NLp5, dashed) for 
the isovector channel are shown. 

[43,44]. 

It has been suggested since some time that the difference in neutron-proton 
flows should be a sensitive observable to study the density behavior of the 
isovector EOS. Several calculations in the non-relativistic framework have in- 
deed shown such sensitivity to the asymmetry-stiffness, but rather extreme 
choices were involved, see refs. [14,45]. In a previous letter [25] we have in- 
vestigated proton-neutron differential flows in the framework of the present 
RMF models. Also in this case a strong sensitivity was found, at high energies, 
clearly due to the Lorentz structure of the isovector interaction. However no 
data exist at present. In this work we will concentrate on pion production and 
spectra. We are also motivated by the fact that recent experimental results on 
isospin effects in pion production are starting to appear, [46,47]. 



4-1 Neutron/proton densities and n /n + ratios 



Observable effects originating from the high density symmetry energy should 
be related to differences on ratios of neutron and proton densities. Thus, we 
consider in the following the ratio of neutrons to protons as a function of time 
and space. One also expects isospin effects on different isospin channels of 
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pions, since they are produced via nn, pp and np collisions through the decay 
of A and N* resonances. 

Fig. 6 shows the time evolution of the ratio of neutrons to protons for regions 
of density higher than saturation, of the multiplicities of A 0, ~- and A +,++ - 
resonances and finally of the ratio of negative to positive pions for central 
(impact parameter b < 2 fm) Au + Au collisions at 0.6 and 1.0 AGeV using 
the NL models without (solid) and with (dashed) the 5 field. In the NLpS 
calculations the neutrons are emitted much earlier than protons from the high 
density phase due to a more repulsive mean field joined to a lower n-effective 
mass. This effect is responsible for the reduction of the n/p-ratio in the residual 
system and, particularly, it influences the particle production, see also [24,48]. 

The four isospin states of the A-resonance are produced from different scatter- 
ing channels, when the threshold energy is available, e.g. nn — > pA~, nA°, pp — 
pA + , nA ++ , ■ ■ ■. Thus A 0,_ (A + ' ++ ) resonances are mainly formed in energetic 
nn- (pp-) collisions. Therefore, the sensitivity on the density dependence of the 
symmetry energy of the n/p-ratio is indirectly related to that of the particle 
production. 

The n/p effective mass splitting in asymmetric matter can also directly affect 
the resonance production around the threshold, with noticeable effects on the 
yields. We remind that resonances will also have isospin dependent in-medium 
effective masses that can be roughly related to the nucleon effective masses just 
using the isopin coupling coefficients in the process A <-> irN [24]. Thus the 
effective masses of the four isospin states of the A-resonance will be different 
when the 5-meson is included in the calculation. In any case we know from 
general grounds, see Section 2, that, in n-rich systems, a isovector scalar meson 
field is leading to a neutron effective mass smaller than proton, in particular 
at high baryon densities. Consequently in the process nn — > pA~ less energy 
will be available for the A~ production in the NLpS case. 

We like to note that this mechanism is not necessarily linked to the density 
behavior of the symmetry energy. E.g., in density dependent coupling models 
we can have at high densities a soft E sym (p) just because the p-meson coupling 
is decreasing but still a large m* — m* splitting if the 5-meson coupling stays 
constant or slightly increasing, like in the DBHF estimations (see Fig.l). 

In the resonance multiplicity (middle figures) one thus sees a decrease of the 
A°'~ isospin states due to the effect of the isovector-scalar 5-meson. The pions, 
on the other hand, are mainly produced from resonance decay according, e.g. 
A~ — > n7r~, A — > pn~ , A + — > mr + , ■ ■ ■. Therefore the decrease of A°'~ reso- 
nances reduces the production of negative charged pions which consequently 
decreases the 7r~/7r + -ratio (right part in Fig. 6). This ratio then appears to be 
sensitive to the isospin term of equation of state at high densities. However it 
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Fig. 7. Same as in Fig. 6, but with the DDHp (solid) and DDH3p5 (dashed) 
models where the isovector couplings are explicitely density dependent. Note the 
larger y— scales compared to the NL cases of Fig. 6. 



is also affected by secondary pion absorption channels which start to influence 
the results at high energies. 

One should note that both NL models used here exhibit a similar asy-stiff be- 
havior at high densities (see Fig. 2). We can deduce that the effective n/p mass 
splitting mechanism described before is very important in inducing isospin ef- 
fects on the particle production [49]. 

Next Fig. 7 shows the same quantities as in Fig. 6 but now using the DDH 
models where all coupling functions are density dependent. In this case we have 
compared the DDHp, very soft symmetry energy, with DDH3p5, relatively 
stiff one. Due to the strong density decrease of the p— coupling function in 
the DDHp parametrization, now the influence of the introduction of the iso- 
vector 5 channel is more pronounced comparing to the NL cases, because of 
the larger differences in the high density symmetry energy and in the n/p 
effective mass splitting. In the transport calculations with the 5 meson (left- 
center panels of Fig. 7) the n/p-ratio and the A~-yield decrease at high density 
is indeed stronger than in the NL models (Fig. 6). Consequently we also see a 
stronger influence on pion production. 

By comparing the Figs. 6 and 7 we see clear differences in the total multiplicity 
of resonances and therefore in the absolute pion yields. One should note that 
the DDH parametrizations differ from the NL models also in the iso-scalar 
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part of the mean field. In particular, the EOS for symmetric nuclear matter, 
which is similar for both models around saturation, has a smaller incompress- 
ibility at p in the DDH cases (see the discussion in Sect. 2) and it is also 
softer at high densities [19,23]. Thus, in the dynamic calculations with the 
DDH parametrizations one naturally obtains more compression in the cen- 
tral region of the reaction which, in turn, is responsible for the production 
of higher multiplicities of resonances. The factor 5 — 7 difference in the A 
formation probability during the compression phase, see the middle panels in 
Figs. 6- 7, seem to indicate that pion absolute yields can be a good probe of 
the nuclear EOS at high densities, apart from corrections due to in-medium 
modifications of inelastic cross sections, resonance decays and 7r-reabsorption 
effects, see ref. [43]. 

The comparison between the DDH and NL models with respect to isospin 
effects depends on several factors. However, the reason of the presentation of 
these different models was to focus on the effects originating from the intro- 
duction of the scalar iso- vector 5 field in the symmetry energy. We have shown 
that there are common features in all models, namely a larger repulsion and a 
smaller effective mass for neutrons in high baryon and isospin density regions, 
for n-rich systems. A decrease of the n/p ratios and of the A~'° multiplicities 
will finally reduce the charged pion n~ /n + ratios. 



4-2 Isospin diffusion 

Before we continue to discuss the pion yields in more detail, we briefly make 
some further comments on the n/p ratios in different density regions. This 
is related to the isospin diffusion from high to low density regions and with 
the subsequent formation of fragments. This has been suggested in various 
works [13,27] as another means of investigating the density dependence of 
the isovector EOS. In Fig. 8 we show the time evolution of the n/p ratio in 
regions with density below (upper panels, "gas and light clusters" ) and above 
(lower panels, "liquid") saturation density for two incident energies. In the 
upper panels we also show the total central baryon density as an indicator of 
the global evolution of the collision. We use the extreme density dependent 
coupling models as in Fig. 7, namely the very asy-soft DDHp (solid) and the 
asy-stiff DDH3p5 (dashed) interactions. 

One sees in the upper panels that there is a strong emission of neutrons to the 
gas just before and at maximum compression for the asy-stiff model, where 
there is a larger differential driving force on the neutrons. The effect is essen- 
tially non existing for the asy-soft model. In correspondence the n/p ratio of 
the liquid (lower panels) decreases below the initial value, more so for lower 
energies. Consequently the chemical composition of the light clusters rapidly 
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Fig. 8. Time evolution of the n/p-ratio (left scale) for density regions below (upper 
panels) and above (lower panels) saturation density in Au + Au collisions at energies 
of 0.25 (left panels) and 0.8 AGeV (right panels). The upper panels also show the 
evolution of the central total densities in units of p sa t (thin lines, right hand scale). 
The models shown are the same as in Fig. 7: DDH — p (solid lines) and DDHSpS 
(dashed). 

formed during the expansion phase will show signatures of the asy-stiffness 
of the EOS, as already suggested in ref.[27]. However, at later times the n/p 
ratio of the dilute region is very similar for the two interactions (slightly below 
the initial value), meaning that at these energies the bulk of the matter finally 
ends up in free nucleons and light clusters. 



4-3 Energy dependence of n /n + -ratios 

We now take up a more detailed discussion of isospin effects on pion yields. It 
was see in Figs. 6, 7 that the 7r~/7r + -ratios stabilize after the compression phase 
and that they depend on the isovector EOS and, particularly, on the inclusion 
of the (5-field. The latter is generally leading to a more asy — stif f EOS 
and in any case to a reduction of the neutron effective mass vs. the proton 
one. The larger neutron repulsion produces less neutron rich matter, which in 
turn produces less negative charged A-resonances. The same effect is directly 
reached by the reduction of the neutron effective mass. Correspondly less 
negative pions are formed, leading to a reduction of the 7r~//T + -ratio of the 
order of 5 to 10%. 
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Fig. 9. Energy dependence of the (7r~/-7r + )-ratio in central Au + Au collisions. Left 
panel: NL results, NLp (solid) and NLpS (dashed). Right panel: DDH results, 
DDHp (solid) and DDH3p5 (dashed). The grey bars correspond to some very 
preliminary data from the FOP I — GSI collaboration [46] . 

All the isospin effects on pion production are generally decreasing with incident 
energy, see Fig.9, due to secondary collisions, i.e. pion absorption and A- 
rescattering. Thus all the isospin dependence of pion production will be on 
average moderated. Such effect was already found in earlier studies using a 
non-relativistic Quantum- Molecular Dynamics (QMD) approach [48]. 

The beam energy behavior of the 7r~/7r + -ratio evaluated at the freeze — out 
time for Au + Au central reactions are shown in Fig.9, for the NL models (left 
panel) and for the DDH models (right panel). We like to note that in both 
cases the 5 influence, difference between solid and dashed curves, is decreasing 
with the beam energy. This is a nice indication of the n/p effective mass 
splitting mechanism in A - production, expected to be more important at the 
threshold. 



For the same systems new very preliminary data on pion production have 
recently been reported from the FOPJ-collaboration at GSI [46], shown as 
grey bars in the Fig.9. These experimental indications are used to test our 
miscroscopic evaluations. 



We see that in our simulations we generally have a qualitative good description 
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Fig. 10. 7r -rapidity distributions for central Ru + Ru collisions with NLp and 
NLp5 models, without the Coulomb field for pion propagation. Full symbols are for 
7r~, open ones for ir + . Data are taken from Ref. [50]. 

of the pion ratio with respect to beam energy, except perhaps at the lower 
energies. The inclusion of a 5 meson in the iso-vector part of the equation 
of state seems to improve the comparison, at least at higher energies. In any 
case it appears that the 7r~/7r + -ratios might be a probe to the role of the high 
density symmetry energy and the role of the virtual 5-meson. It would be nice 
to look at more exclusive data on pion production in asymmetric systems. 
This will be the subject of the next section on pion flows. 



5 Isospin effects on pion flows 

We now discuss in more detail the velocity distributions of pions, i.e. flow 
results. In spite of the relatively low asymmetry of the system (N/Z = 1.18) we 
will study central 96 Ru+ 96 Ru collisions at about 1.5 AGeV, since recent data 
are available from the FOP I collaboration [50] . We will show symmetry energy 
effects and how in this respect the Coulomb interaction in pion propagation is 
important. In particular we can expect the role of the Coulomb field to be not 
negligible in high density regions, for the interaction between the produced 
pions and the charged hadronic enviroment. We will also see that without 
taking the Coulomb effects into account, the differences originating from the 
S meson are largely reduced. 
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Fig. 11. vr^-rapidity distributions for the same system as in Fig. 10, for the NLp 
(full circles) and NLp5 (full squares) models, but including the Coulomb field in 
the pion propagation. Upper panel: tt~ . Lower panel: ir + . Data (open diamonds) 
are taken from Ref. [50]. 

In Figs. 10 and 11 we show the rapidity distributions for in comparison to 
the data for the NLp and NLp5 models. The normalized rapidity is defined 

by Y = {Y/Yproj)cm- 

In the calculations of Fig. 10 the Coulomb interaction was turned off for the 
pion propagation. It is seen that without the Coulomb interaction the rapidity 
distribution is generally too high compared to experiment, particularly so at 
mid-rapidity. Including the Coulomb, Fig.ll, the distributions are reduced, in 
much better agreement with the data. The effect of including 5-fields is also 
enhanced. 

During the high compression stage the Coulomb field accelerates the positively 
charged pions away from mid-rapidity. The opposite is happening with the 
7r~'s that will be attracted towards the higher density region where they will 
suffer a larger absorption. This leads to a broadening of the pion distributions 
and to a general reduction due to the enhanced rescattering and absorption. 
Similar effects were also seen in the non-relativistic calculations of ref. [48]. 

The combined role of the 5- and Coulomb- fields on the difference between 
positively and negatively charged pions is a more subtle effect. As we have seen 
previously the structure of the symmetry term affects the isospin diffusion from 
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high to low density regions. Moreover the S contribution is directly reducing 
the negative A-resonance production through the n/p-mass splitting. As a 
consequence the NLp5 model leads to a more positively charged hadronic 
central region and to larger Coulomb effects, The freeze-out it~/it + ratio at mid 
rapidity is then reduced. However, the Coulomb field, acting more repulsive, is 
inducing larger reabsorption effects for high momentum tt + . As a consequence 
we can have an enhancement of the 7r~ /tt + ratio for high-rapidity pions. The 
influence of the 5-channel is thus observable in terms of pion stopping. 

This can be seen in Fig. 12 (left panel), where we show the rapidity dependence 
of the 7r _ /7r + -ratio for the two NL models. For mid-rapidity the inclusion of 
the 5-field reduces the ratio with the mechanisms described before. Since the 
mid-rapidity pions are the most abundant (Fig. 11), this is consistent with the 
global reduction of the 7r~/7r + -ratio (see Fig. 9, left panel). 

5.1 Transverse momentum dependence 

All the effects discussed before should be particularly important in the com- 
pressional phase of the collision, while pions are produced all over the reaction 
evolution up to the freeze-out. A good selection procedure will be to look at 
particles with high transverse momenta P t . It is indeed physically intuitive, 
and confirmed by simulations [25,48,51], that particles emitted from the earlier 
high density phase have higher P t 's. 

The P t - dependence of the 7r _ /7r + -ratio in central Ru + Ru collisions for mid- 
rapidity pions is shown in the right panel of Fig. 12. It is seen that the inclusion 
of the 5-field strongly enhances the ratio for high-Pi pions, bringing it into more 
agreement with the data. 

It is here that one sees the combined symmetry and Coulomb effects most 
clearly. With inclusion of the 5-field the high density phase is relatively more 
rich in protons and positively charged resonances. The Coulomb force accel- 
erates 7r + away from this region and so reduces their number at mid-rapidity. 
All that then enhances the 7r~/7r + -ratio for high-Pi particles. 

We have repeated the analysis for the two DDH interactions already studied 
before, DDHp and DDH3pS. The inclusion of the ^-channel leads to the same 
effect. We see an increase of the 7r~/7r + -ratio for high-Pi pions at mid-rapidity, 
as observed in experiments. 

We have shown that 7r~//T + -ratios for high transverse momentum pions are 
a sensitive probe to the structure of the symmetry term, in particular to the 
contribution of a 5-field. The trend of the data with rapidity and transverse 
momentum is affected by the dynamical effects of the inclusion of a isovector 
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Fig. 12. 7r~/7r + -ratio for central Ru + -Rit-collisions as a function of rapidity (left) 
and transverse momentum (right) at mid-rapidity (|l" ( - - ) | < 0.15). Same models as 
in Fig. 10. Data are taken from Ref. [50]. 

scalar field in the hadronic model. This seems in better agreement with the 
data, in particular for high-P t emitted pions. 

We finally remark that this already clean indication could be further supported 
by studying instead of m Ru + 96 Ru (N/Z = 1.18) more asymmetric cases, like 
132 Sn+ 124 Sn (N/Z = 1.56) or 197 Au+ 197 Au (N/Z = 1.49) colliding systems, 
and looking at similar more exclusive data. 



6 Summary and outlook 

In this paper we investigate in detail the density and momentum dependence 
of the symmetry energy in nuclear matter properties, finite nuclei and, in 
particular, in dynamical situations of intermediate energy heavy ion collisions. 
We discuss isospin effects in terms of the Lorentz structure of the isovector 
mean field. In addition to the usual vector-isovector p-meson we introduce 
a scalar, isovector 5-like field, naturally keeping the same bulk asymmetry 
parameter at saturation. This will, however, change the high density and high 
momentum behavior of the symmetry term. We study the effects arising from 
the 5-meson using different models of the Relativistic Mean Field (RMF) 
theory, then stressing the role of Lorentz covariant properties of the isovector 
mean field in a general way. 
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The comparisons between Non-Linear Walecka type (NL, constant couplings) 
and the microscopically motivated density dependent hadron field (DDH, 
density dependent couplings) models show that the effects of the 5 meson 
are of general nature and that they originate from the different covariant 
properties of the vector p— and scalar 5— fields at high baryon densities. 

We have applied the models of asymmetric nuclear matter to finite nuclei and 
heavy ion collisions. 

The studies of nuclear structure show that the effects of including the isovector 
^-channel are very small, except if one goes much beyond the valley of stability 
near the neutron drip line. Observable effects are found for very neutron-rich 
nuclei as an increased neutron diffusion at the nucleus surface and a larger 
binding energy per nucleon. Both results are due to the fact that in n— rich 
systems at low densities neutrons feel a larger attractive scalar self energy when 
including the 5 meson. However, these effects are only moderate, in agreement 
with similar analyses of ref.[23], making it difficult to draw definite conclusions 
on the effective interactions. The main reason relies on the observation that 
in finite nuclei one explores the EOS around and below saturation density. 

Thus, in order to investigate the Lorentz structure of the high density sym- 
metry energy we have applied the models to dynamical situations in heavy 
ion collisions at intermediate energies, where relatively high baryon density 
regions can be reached during the reaction dynamics. 

The isovector scalar and vector 5- and p-fields cancel each other in a similar 
way as for the isoscalar a- and cj-fields. Thus, with the constraint to have 
the same symmetry energy at saturation, this cancellation is less effective at 
higher densities, with the repulsive p-meson gaining more importance [15,16]. 
With this pure relativistic mechanism the inclusion of a 5-field will make 
the symmetry energy more repulsive at high baryon densities. Another very 
important transport effect of a scalar field in the isovector channel derives 
from the influence on the effective masses. In asymmetric matter we have now 
a splitting of the n/p effective masses, also increasing with baryon density. This 
will affect nucleon velocity distributions and particle productions, particularly 
around the threshold. Finally, at high relativistic energies, due to the large 
velocity fields, one can expect large effects directly from the different covariant 
structure of the p and 5 mesons [25]. 

We have studied heavy ion collisions in a relativistic transport model of Boltz- 
mann type and analyzed the results in terms of neutron and proton densities 
in space and time, particle production and pion flow observables. Within the 
same frame we have studied nucleon flows in a previous letter [25]. 

For n-rich colliding systems we find a strong decrease of the n/p-ratio in the 
bulk matter when including the 5 meson in the isovector mean field. We in- 
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terpret this effect by the joint action of a stiffer symmetry term and of a 
decrease of the neutron effective mass, more pronounced at high baryon den- 
sities, which both lead to more repulsion for neutrons than for protons. The 
influence of the 5-meson in the dynamical situations has been mainly studied 
within the iVL-models, where the differences in high density symmetry energy 
are actually not so strong. However, the DDH models of density dependent 
coupling functions show very similar effects, even more pronounced due to the 
softer character of the symmetry energy. Moreover the dynamical implications 
of nucleon effective mass splitting are present in all cases. Therefore, the in- 
fluence of the <5-channel is a general feature arising from the differences in the 
Lorentz structure of the high density symmetry energy. 

The relativistic form of the high density isovector EOS can also be directly 
seen in terms of particle production which can be experimentally measured. 
For this reason we have focussed our study on the influence of the 5-meson 
on pion production. We evaluate from the reaction simulations the ratios of 
charged pions as function of time and space, energy, rapidity and transverse 
momentum. In general for all the colliding systems and energies we observe a 
reduction of the (7r~/7r + )-ratio in the models containing the 5-meson. 

This is due to the fact that the 5 meson makes the high density region more 
proton-rich, see before. This induces a reduction in the formation of A ' - 
resonances and thus in the production of negative charged pions. Moreover 
the lower neutron effective masses represent an efficient mechanism to directly 
quench the formation of negative A-resonances. This is particularly evident 
at lower energies, around the threshold for resonance production. 

The effect of the 5-field as discussed above is enhanced in the momentum 
spectra, i.e. the rapidity and transverse momentum distributions, due to the 
interplay with the Coulomb interaction which affects differently the n + and 
7r~ distributions. Thus a clear, and likely detectable, effect of the 5-field is pre- 
dicted in the flow observables of the (7r~/7r + )-ratios, as well as in the incident 
energy dependence. 

The comparison with the experiments seems to indicate that the introduction 
of the S meson in the isovector mean field is important, i.e. as a way to promote, 
in n— rich systems, a larger vector repulsive field for high energy neutrons, joint 
to a smaller effective mass. We remark that both effects are mainly due to the 
relativistic dynamics and not necessarily directly related to the high density 
behavior of the symmetry energy. 

However, in order to strengthen these conclusions one needs more precise ex- 
perimental data, with more asymmetric, even unstable, ion beams, not only on 
pion production but simultanously on collective flow of neutrons and protons 
to study the so called isospin collective flow [25]. 



26 



In the same heavy ion collisions, with high isospin density, the threshold pro- 
duction of other particles and resonances would also provide very relevant 
information. 

At this level of investigation we conclude that relativistic studies of asymmetric 
static and dynamic nuclear matter appear very important in order to directly 
probe the Lorentz structure of the isovector part of the nuclear equation of 
state. 
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